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‘The Aerospace Synthesized Aerodynamic Heating Program is written in the FortraniV 
language for the IBM 360-65 computer. This program calculates heating and internal 
conduction for shapes approximating those of missiles. To do this, it couples various 
correlations for external high speed convective aerodynamic, for nonaerodynamic, heating 
with internal transient heat conduction and surface ablation. These may be applied-to flat 
plates, wedges, cones, sphere-cones, spheres, or cylinders with laminar and/or turbulent 
heat transfer equations applied to the configurations listed. The method is sufficiently 
general to allow many variations in the physical properties of materials and boundary 
conditions that one might wish to consider. ( ~~ 
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Details of the method of analysis are covered in Section 2, while a program listing is — 
in Appendix B, along with sample data cards. 


h 


Sc AN RSA Ti, AMMEN at ANC GA RMN iN SSO ALA AI el MURR ie sods Lk sense AEA ot gel pel italia es Walid TM tei a 


a tu, 


ote 


DD r..1473 UNCLASSIFIED 


i: Ain ohari aba 


( 
I 
i 
i 
| 


| 


ESET EANNE ENP AY NN ATTN NNT NET EN EEN TT NBN PAT A EIA NT TB HT EST ET oP TL 7 


4 i 
APSO Be ama a a 
Te Sen eget me ” cme NCR MLM HR anmeoenn | em 


Won 
ANIA ahi nn ae oro 
NC ORO Roe AEN A hana ra 


ROLE 
UNCLASSIFIED 
Security Classification 


assification 


ty 


high speed convective heating 


aerodynamic heating 
multi configurations 


: 
: 
; 
Z 


surface ablation 
heat conduction 
eolar heating 


HPT RO RR 


aged 


Pk SOU EIT] Fa bt eu 


aN Pd 
wy 


ert 


[Ss 


A ge a RR 


Air Force Report No, Aerospace Report No. 
SAMSO-TR~ 72-181 TR-0172(S2816-75)-1 


AEROSPACE SYNTHESIZED AERODYNAMIC HEATING PROGRAM 


Prepared by 
M, J. Gyetvay and J, M. Kohlenberger 


72 MAR 15 


San Bernardino Operations 
THE AEROSPACE CORPORATION 


Prepared for 


SPACE AND MISSILE SYSTEMS ORGANIZATION 
AIR FORCE SYSTEMS COMMAND 
Air Force Unit Post Office 
Los Angeles, California 90045 


Approved for public release; distribution unlimited, 


= Se ee 


° 


es 
23 
34 
ia 
.— 
zs 
P= 
ff 
= 


Af 


cen teammreis + 4 einen ett ete RA RE OO St NLU TCL ANIC CETTE et ARE Ae RE TAD Tt, 


FOREWORD 


This report by The Aerospace Corporation, San Bernardino Operations, was 
prepared under Contract No, F04701-71-C-0172 as TR-0172(S2816-75)-1, The 
Air Force program monitor is Capt. T. W. Swartz. The dates of research for the 
report include the period July 1963 throvvh April 1971, This report was submitted 
by the authors April 1971, 


The method presented in this report for solving in-depth material heating was 
first proposed by D, J. Smith and L, Fisher (Ref, 1) of Convair. A preliminary 
program was formulated at Aerospace by R, Chambers, M, Kausch, and 
J. Kohlenberger in 1963. The authors have since developed the program into its 
present form, based on the work of the above named people, 
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AEROSPACE SYNTHESIZED AERODYNAMIC TR-0172(S2816=-75)-1 
HEATING PROGRAM, by M. J. Gyetvay and 72 MAR 15 
J. M. Kohlenberger 


PRT ene Rom nON 


The Aerospace Syrthesized Aerodynamic Heating Program is written in the Fortran IV 
language for the IBM 360/65 computer. This program calculates heating and internal 
conduction for shapes approximating those of missiles. To do this, it couples various 
correlations for external high speed convective aerodynamic, nonaerodynamic heating, 
with internal transient heat conduction and surface ablation. The heating may be com~ 
puted for flat plates, wedges, cones, sphere~cones, spheres, or cylinders with laminar 
and/or turbulent heat transfer equations applied to the configurations listed. The 
method is sufficiently general to allow many variations in the physical properties of 
materials and boundary conditions that one might wish to consider. 
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Details of the method of analysis are covered in Section 2, while a program listing is 
in Appendix B, along with sample data cards, (Unclassified Report) 
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NOMENCLATURE 

a 
A’ = conductance At = time step 2 
Cf = friction coefficient a = angle of attack 3 
h = heat transfer coefficient € = emissivity F 
a = coefficient o? pressure 6 = cone 1/2 angle 3 
ae = heatiag multiplication daction A = sweep angle g 
Us = entiasy ratio p ~—-=_- viscosity F 
BK = eoedh -ivity p = density 


= iat ecefficient based on enthalpy SUB AND SUPER SCRIPTS 
= enthalpy 
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abl = ablation 

b = base temperature (400 °R) 

eff = effective 

i = calculated from reference enthalpy 


n = a particular nodal point, or time 
step 


H 

I 

J = conversion facto 
M = Mach number 

Pp = pressure 

Pr = Prandtl number 
Q = heat transfer 

Q = heat transfer rate 7 
R 

Re 

St 

T 

U 

x' 
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r = recovery property 

= radius of nose 8 = stagnation property 

= Reynolds number sl = sea level 

= Stanton number sph = sphere property 

= temperature tan = tangency poi: 

= velocity Vv = vaporization 

= diameter of sphere w = _wall property 

= lor nesiertag heating factor as =, = edge of boundury layer 
vA = fictitious wetted distance : Shinavaas 

(see eq. 18) 0 = free stream property 

a = nodal thickness 
cp = specific heat 
g = gravity (32.2) 
r = recovery factor 


= wetted distance 
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SECTION I 
INTRODUCTION 


T..2 Aerospace Synthesized Aerodynamic Heating (ASAH) program is written 
in the FORTRAN IV language for the IBM 360/65 Computer, itis a flexible computer 
program that couples various correlations for external high speed convective aero- 
dynamic heating with internal transient heat conduction, The flexibility of the program 
is derived from the many modes of operation available to a user and the ability to run 
parametric studies in a Simple manner, A programmer may select from the following 
list of body configurations: flat plate, wedge, cone, sphere, cylinder, and nose blunting, 
Ap option is available for either ablating or nonablating surfaces, Laminar and/or 
turbulent heat transfer equations are available for each of the above configurations, 

Boundary layer transition is based upon either an input transition Reynolds number, an 
input transition altitude, or in the absence of an input altitude or Reynolds number 
, value, the program will compute its own transition Reynolds number baeed upon a 
0) correlation by McCauley (Ref. 2). 


For input, the program requires specification of a trajectory, physical and 
‘thermal Properties, and control instructions, The trajectory input spec ‘fies altitude, 
velocity and angle of attack as a function of time (Fig, 1~1), Physical dimensions such 
as distance from forward edge (Fig. 1-2) and material segment thickness, and thermal 
properties such as material specific heat, density and conductivity must also be 
supplied, Specific heat, density, and conductivity may be inserted as functions of 
temperature, Control instructions define modes of cperation such as configuration 
* identification, transition criteria, boundary conditions, and computational time interval, 
- The program will compute an optimum integration interval, after the first iteration, 
All inputs are entered by way of an IBM input routine called NAMELIST, which allows 
a large degree of freedom as to input format, Heating can be either internally computed 
: e aerodynamic, solar, or input us a function of time, Inputs are stored and kept until 
changed, thus facilitating parametric studies, 
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The printed output for this program includes edge of boundary layer properties, 
reference conditions, free stream conditions, heat flux, dynamic pressure, Reynolds 
number, shear, material ablated when applicable, and temperature distribution through 
the heat protection system, When the nose blunting option is used, radius of a blunt 
nose, stagnation recession, sonic point heat rate, shear, and heat of ablation for a 
two-component ablative composite is also provided. 
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SECTION II 
COMPUTATIONAL PROCEDURES 


The actual flow of computation through the ASAH program is dependent. on 
the particular problem being solved. A generalized discussion of the flow of typical 
on 
computations and the subroutines involved will be given in the following paragraphs. 
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The variable, time, is used as the independent variable to control the overall 
flow of the program, Subroutine TRAJ linearly interpolates between appropriate 
trajectory inputs to corresponding free stream velocities and altitudes. These inputs 
are used in conjunction with subroutine ATMOS to provide the free stream atmospheric 
properties for any altitude up to 295, 278 feet (Ref. 3). ‘ 

When the wedge or cone shape configurations have been specified, the program 
algebraically adds the angle of attack to the body half-angle to determine if the SHOCK 
subroutine or PME (Prandtl-Meyer-Expansion) routine should be used. The SHOCK 
subroutine uses a similarity parameter to determine the edge of boundary layer stream 
properties behind the shock. The PMF subroutine utilizes the Prandtl-Meyer isentropic 

relationships to determine fluid properties aft of a sharp corner. 


Next, the program computes density, viscosity and Reynolds number. The 
program uses Eckert's reference enthalpy to evaluate high Mach number heat transfer, 
These calculations are basically carried out in, or under the control of, subroutine 
THERMO, The Aerodynamic heating calculation for flat plates, wedges and cones arc 
carried out in subroutine ECKERT. Subroutine ECKERT uses the modified Reynolds 
analogy to correlate skin friction with heat transfer. In the case of laminar flows, 
skin friction is evaluated by the Blasius relationship. When the flow is turbulent, skin 
friction is determined using the Schultz-Grunow relationship. The aerodynamic 
heating calculations for spheres, cylinders, and blunted cones are handled in 
subroutine KEMP, using the Kemp-Riddell stagnation correlation (Ref. 4), or the 

Sibulkin equation for Mach numbers less than 3. There are also options for calculating 
heating inputs using the Hankey sonic point heating formulation, in subroutine HANKEY. 
Another option is available in subroutine CONSTQ. This routine provides heating values 


based on interpolation from tabular input. 
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Upon completion of the external heating calculations, the program calls 
on subroutine ENERGY, This routine performs an energy balance computation at the 
flow field ~ materia! interface, 


The thermal protection system (heat sink or ablating surface) is divided into 
a number of segments, With the net heat load as a boundary condition, Fourier's 
transient conduction heat transfer equation is solved in finite-difference form. The 
calculations are carried out under the control of subroutines ONED or TWOD for either 
one or two spatial dimensions, 
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The one~dimensiona! option allows the heat protection system to be divided into 
a maximum of 49 segments. The two-dimensional option allows the system to be 
characterized by a 20 x 20 matrix of points, Material specific heat and conductivity 
can be inserted as functions of temperature, Each matrix element or segment may 4 
have different material properties, The program obtains particular values of the 3 
material properties, for use in computation, by accessing subroutine GETCPK or 
GPCPK2, These routines access the proper input material properties tables as a 
function of temperature and provide the needed values, } 
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When the ablation option is exercised, the energy balance at the body surface 
allows for energy to be absorbed by the ablating material, This adjustment is carried 
out in subroutine ABLATE, A gross effective heat of ablation is computed as a function 
of heat of vaporization and enthalpy difference, The rate at which material is removed 
q from the surface is determined by the difference between net heat flux to the surface 
q and the heat conducted away into the next material segment, When the surface 
temperature reaches the ablation temperature and net heat flux to the surface is 
positive, the ablation is performed, 
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The nose-blunting option uses the Kemp stagnation point heating and engineering 
correlations to obtain heating at the sonic or tangency positions. These calculations are 
carried out under the control of subroutine ABLUNT, 
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The computations that have been, so far, described comprise those required to 
advance one step in time. The main program now performs various checks relating 
to the print and punch options. It then does a stability check to select a proper 
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3 3 time step fox the next step in the computation, Control is then transferred to the 
7 beginning of the computational loop and reenters subroutine TRAJ. This mode of 
operation is continued until such time as the case is completed, or an unrecoverable 
problem is encountered, When this occurs, the program transfers to a call to sub- 
routine INPUT to look for subsequent cases, 
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SECTION IN 


METHOD OF ANALYSIS 


This section of the report reviews the various methodologies that are used by 
the computer program to handle the individual aspects of the overall! prob!'em, The 
intent is to document the equations and data that are in the program and indicate the 
source of major portions of the formulation, No attempt is made to derive the equations 
used or to justify the assumptions inherent in the computer program, 


3.1 HEATING RATES FOR FLAT PLATES, WEDGES AND CONES 


The reference enthalpy method proposed by Eckert (Ref, 5) is used to compute 
the convective aerodynamic heating rates, This method is based on the assumption that 
the incompressible mass, momentum, and energy equations can be used for compressible 
flow solutions, providing the transport properties of the gas are known and evaluated at 
a reference enthalpy. Real gas effects, including gas disassociation and ionization, 
must be taken into account also in the determination of the thermodynamic properties of 
the gas just outside the boundary layer, 


3.1.1 Laminar Boundary Layer 


The Blasius solution for an incompressible boundary layer is 


St = a pr 2/3 (1) 
and 
cr= 02864 (2) 


(Re 


where St is the Stanton number, while Pr is the Prandtl number (= 0.65) and 
Re = pUX'/y is the Reynolds number. Using Eckert's reference enthalpy 


* t 
Re* = os 4 (2a) 
Hi = 3600 St P35 Us 4 (3) 
3-1 
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where p and yu are evaluated at a reference enthalpy, I*, and U is evaluated 
just outside the boundary layer, and g is gravity. 


When computing recovery enthalpy the recovery factor is Pr" with n= 1/2 
for laminar flow and n= 1/3 for turbulent flow. Since Pr was assumed = 0.65 
the following values for recovery factor (r) are used: 


0.866 for turbulent flow 
0. 806 for laminar flow in nonstagnation areas 


1,0 in stagnation areas (sphere and cylinder). 
Q = Hy dh) 


it 


‘= 


(4) 


The thermodynamic properties of the air ahead of the shock, Ps Ty 


and the free stream velocity and Mach number, UL, and M,.» can be determined 


from the given flight conditions. An approximation of the ARDC 1959 atmosphere is 
used, The free stream and wall enthalpies can be determined once the ambient 
temperature and wall surface temperature are known, The thermodynamic properties 
just outside the boundary layer are determined from the properties given by Schadt 


and Livett (Ref. 6) once the hypersonic similarity parameters, M,, Sina and M.. Sing 
are known (a, here is total of a +6). Eckert states that 
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= 1, +0.5 (I, - 1) + 0.22 (I, - 5) (5) 
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where 
ru; 
I. -L+3y z (6) 


‘.ne transport properties are given by Hansen (Ref. 7). Curve fit equations were 
written making p* and p* dependent on temperature. The terms Py and M, are 
evaluated at a base temperature equal to 400°R suggested by M. F. Romig (Ref, 8). 
For cones the following correlations are used to define edge flow field properties: 
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when M.. Sina < 4,0: 


(us es Us) x10 = -0,107 (M,, Sin a)® + 1.91 (M Bin ay” 
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+ 0.04 M, Sina + 1.57 (7) 
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when M,. Sin a> 4.0: 


(we U?) x 107° = - 0.0062 (M,, Sin a)® + 1.414 (M Sin a)” 


+ 1,038 M. Sin a + 6.187 (8) 


with M Pia U,,/1086. 


Pp I 
The following hypersonic similarity equations for = and 2 are then employed: 
ie) [oe] 


when 1.5<M, Sina < 5.0 


ere M_Sina)® +1.06 (M_ Sin a)” 
Bp” * ( a Qa 96 ( oo Hn a) 
+1.16 (M,, Sin a) (2a) 
when M,, Sin a > 5.0 
BO sa aargiae Sin a) + 0.8373 (M_ Sin a)” 
P- - ( i na) ° ( co Sin a) 
+4.177 (M,, Sin a) - 3.182 (9b) 


and when M, Sina < 8. 


2 
I; = Ty {1.03 + 0.0827 My, Sin a + 0.2174 (M,, Sin a) 


3 (10a) 
~- 0.0006956 (M. Sin a)} 
or, when M, Sina> 8. 
a . 2 
I, = 1, [1.106 - 0.3685 M, Sin a + 0.3466 (M, Sin a) 
ho C. 
007766 (M,, Sin a)*| Aa 


with M Be U,,/1036. 


Similar correlations are used to compute edge of boundary layer fluid properties 
for a wedge geometry. In the case of the flat plate PS and U,, are used in place of 
Ps and U 63 the ambient air properties are assumed at the edge of the boundary layer. 
The density is obtained from the ideal gas equation, and the viscosity from the Sutherland 
equation. The heat transfer parameter obtained by Eq. (3) is multiplied by a factor of 
1.73 for conical flow with laminar conditions (Mangler transformation). 
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3.1.2 Turbulent Boundary Layer 
The Schultz-Grunow correlation for an incompressible turbulent boundary 
layer is used with the modified Reynolds analogy to predict turbulent heat rate at 


high velocities 


st= GE p,-2/8 (11) 
with 
Cf = 0.370 55 (12) 
{log (Re*)/2.3031“° 
The equation for the heat transfer parameter can be reduced to the form 
H, = 3600 g St Us Ps (13) 
and 
Q = H(,-) (14) 
p* U, x' 
Re* = oo (15) 


For a cone the heat transfer parameter is increased by a factor of 1.176 (Mangler 


transformation) over the flat plate heat transfer rate. 


Hypersonic similarity parameter correlations yield I 6? Us ’ Ps - Other 
fluid properties at the edge of the boundary layer, are evaluated with the following 


real gas correlations: 


when I5< 1300 


p5= 5.76 x10 1, 9-849 Pp (16a) 
yg = 4.2642 x 107 1, 9° 498 (6b) 
when =I, 2 1300 
: 6. ~0,584 
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7, 0,228 


= 2,8428 x 10 5 


I 


Ms 


Re, = Usp5X'/its | (18) 
This locul Reynolds number is tested against a transition criteria to determine 
laminar or turbulent flow. 


For laminar and turbulent flow, X' the input wetted-distance, is used to 
compute local Reynolds number. 


3.3 SPHERE-CONE 
8.2.1 Laminar Flow 


When the flow over a sphere-cone is laminar, the heating correlation is 
identicelly that recommended by Griffith-Lewis (Ref. 9). 


Local flow properties are evaluated at a nondimensionalized length: 


2 
_9 Cos@x 
a= "0.802 R (19) 
The following array of algebraic equations are solved to determine a local 
heating rate. The factor Y is determined from table lookup and is identical to the 
valve, suggested by reference 9. 


I = 0,24T (20) 
[9] co 
I I 
6 ° 
BeS (+s 22) (21) 
2.581, 0.253 
Cs2 = 0.588 (22) 
w= 0.173 x 107 + 0,00833 x 107° I (23) 
[+ e) foo) 
64.4 Po UL R 
OT. ee ere (24) 


& ry 


Then, if Z< 1.5 the coefficient of pressure is defined as 


2 


Cp = 2.4660 (25a) 


or, if Z >1.5 


Cp = 267 (-6.32° + 29,327 - 52,2 Z° + 43.152" - 15.52 + 2.76) (25b) 
MYO" Cs2 (1.0+ ee 
M Cp 
Chet ee ae (26) 
0,632 / Re 
Surface heating is then computed as 
1,2288x10° ChP M2 
H, = a 2 (27) 


wo 


Finally, the convective heat transfer, to a surface of unit area, is computed by the 
equation 


Q =H, (I, - 1) (28) 
3.2.2 Turbulent Flow 


For turbulent flow on a sphere~-cone, a Schultz-Grunow heating solution is 
employed. The flow field properties are computed by assuming isentropic expansion 
from the stagnation point to a given pressure. This pressure is obtained from the 
correlation suggested by Griffith and Lewis. The local flow field properties are 
computed as follows: 


: 2 
Po=P(0.7M_Cp+1.0) (29) 


0. 286 
0.8333 ue, 1,0 - =) 1, 168 M2 - 0. 1668 


(30) 
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at this point Eckert's reference enthalpy, Eqs. (5) and (6) are used along with Eqs. (11) 
a through (17). 


The turbulent heat transfer coefficient to a cone is determined by increasing 


the flat plate correlations by (2yt/ A Finally the convective heat transfer to a surface 
of unit area is computed from Eq, (27). 


3.3 SPHERE AND CYLINDER 

4 The heat transfer to the stagnation point of a sphere given by the Kemp and 
RO Riddell equation is used: 

ROR 3.15 

po EO g = 8284 x10° | % [%s (i, ~ hy) ag 
BOE sph \ 1 \26000 (i, =129.5) 


PR where p,, is sea level density. 


E: : An approximation to the Sibulkin equation (Ref. 10) is used for M,, < 3.0 


e O& [Pog Y 
a Qeph = 174-0 = 2 (i, - i) (34) 


On an unswept cylinder, the convective heat transfer to the stagnation point is 
; Qeyl =Q 'sph x factors (35) 
¥ The following equations are used to account for sweep angle effects 
4 _ 1.5 
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where 


(37) 


A a7 Tan72 cosa _sinA | 


i . V1.0 - (cosw sinA)” 


A program option is available whereby Q; a heating input, can be found by way of 
a table look-up. If Cq is zero or negative, the program looks up the proper Q; 
and sets Cq = 1 at which point QNET is found by the ordinary external energy balance 


used for all situations. 
Ts 4 
0.174€ Too (38) 


3.4 TEMPERATURE DETERMINATION 


Pe ee St 


4 QNET = (Cq Q)) - 
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The material temperature is a function of both the net heat transfer to the ; 
iw wall (i.e., including radiation losses) and the diffusivity pr~perties of the material. 
In order to evaluate the surface temperature and the material temperature gradients, 2) 
the so called "lumped parameter" method is employed. The following section is a 
general description of this method. 


that 


The material is divided up into an arbitrary number of small segments, 

| both perpendicular and parallel to the flow direction as shown in Figs. 3-1 and 3=2. 

Then for some small time increment (controlled by stability considerations) the net heat 

flux to the surface is determined. Knowing this value, the net heat transfer from 

segment (1) to all the adjacent segments is computed for the same small time 

increment, In this calculation it is assumed that all the mass is concentrated in a 

point source at the center of the segment and that the heat transfer takes place over 

an area equal to the interface between the segments and a distance equal to the spacing 

between the centers of the segments. The difference between the heat transferred to : 
segment (1) by convection and radiation, and that which is transferred to all the 
adjacent segments by conduction causes 2 certain tempezature rise in segment (1) 
during the calculation time interval. By making the calculation time interval and the 
size of the segments very small, the accuracy of the method may be made as good as ; 
is desired within the accuracy of the computer being used. a ) 4 
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A forward marching explicit method of computing the temperatures is used 
where the temperature of a given segment at the end of a time interval is dependent 
upon the temperature of the given segment and adjoining segments at the beginning of 
the time interval. 


3.4.1 One-Dimensional Conduction 


The material or materials may be divided into n segments, as shown in 
Fig. 3-1, where n<40, They are numbered from the outer surface inward, 


hen BB gs 1S DRS MW REME AC Seo heathen a ibe ote bd adr ase a ba a 


The heat balance for conduction may be written as follows: 
for internal segments, 


T -T 
T= T +At—och ey 
n n 


1/*n , 2) 
2 K, K-1 


c 


ee a ee ee Lost 


T 
i a BS Se (39) 
1 a rs ane 
2\ K Kael 


for the surface nodal point (n = 1) 


eo ORES Aa tetany MTA ata a + 


Tae Aa (is fen) (20) 
= + 


Where a is the temperature at segment n for the succeeding time increment, "a" 
is the segment thickness, and K is conductivity. The conductance is defined as: 


At 


AMS 
25. 0 Py Py By 


. (41) 


For the last segment a quantity HINSD = 12 (TINSD ~- tT) is used to account for 
internal heating (if present). 


3.4.2  ‘two-Dimensional Conduction 


8 The two-dimensional heat conduction technique involves setting up a structure 
: in a matrix such that the maximum number of rows = I and maximum number of 
columns = K. The material must be homogeneous in the X direction but not 
necessarily in the Y direction. This is shown in Fig. 3-2, ‘ 
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The first slab in successive rows is numbered (1, 1), (2, 1),.. (1,1), there 
are IxK nodes, or elements, therefore the temperature at a node i, k, time t + At, 
is dependent upon the temperatures of the contiguous segments at time t, The general 
picture for a given nodal point is shown in Fig. 3-3, 


71 26125 02 


Figure 3-3, Nodal Subscript Relationships 


As shown in Fig. 3-3, the equations for finding the temperature Th for segment 
k 


(i, k) are: i, 
At 
RPRQOR eee (42) 
5 
25 Ps OP; 4 & By 
b, CC 
= k ~~“i-1,k 7 
TERM 1 = Tak Ti) (43) 


7 fis ea 

2\Kiaik Si-tk 

where K is conductivity, and CC is a multiplying constant that allows configuration 
flexibility, 


(CC = 1 if conduction desired and 0 if no conduction), 
and (Thy 7 Ty, 4) = heat from (i-1,k) segment (44) 
) 9 


TERM 2 in like manner with heat from (i+1,k) segment 


8-12 


C3 
AN Lc Sia Dak Di a al ll lhe, Ll ancl EEO RN Nin et ON mE cans te Blaha eee et ALN OTS 
Dia hye Mas Bad l # ate a UPR RET eT phere Le cat! Bde tl te ptt wy gat eh Ba Mea bi ica kot: desea g! 


ar ia 4 lt NOON MN LD LDA Naish aa NE beset emace eietuemitnte mae wnremcnaem 
aba ly hE Se tn fei Datei tte Da reese! FUMES Shu deSt te ak UR gOS sig te wig ul 


TERM 3 in like manner with heat from (i, k-1) segment 


TERM 4 in like manner with heat from (1,k+1) segment 
(44) 
A = FACTOR(TERM1 + TERM2 + TERM3 + TERM4) 
i,k 


3.5 ABLATION TECHNIQUES 
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Some of the more promising types of heatshield materials for high speed 
flight include those whose surfaces melt or ablate because of the aerodynamic heating. 
In this kind of process it is necessary to compute the rate of surface ablation. This 
is done in the following manner. The heat transfer to the ablating surface is calculated 
as previously until the surface segment reaches ablation temperature. At this time 
two simplifving assumptions are invoked: 


(1) The segment temperature is maintained constant and equal to the 
ablating temperature. This is usually considered to be the melting 
temperature of the fibers (glass, quartz, asbestos, etc.) rather 
than the melting temperature of the binder. 


wecwtar 


(2) The integrated heat flux to the segment is compared against the 
effective heat of ablation of the material, The effective heat of 
ablation is determined from test data and numerically represented 
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by the following equation: 
Hap = Hy * 7 (E - AY) (45) 


Hy and 7 are input parameters to the code. 


Material weight loss is determined by integrating in time. 
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When weight loss equals equivalent weight of a segment 
thickness that nodal point is removed from the solution 
temperature matrix. 


3.6 NOSE BLUNTINC TECHNIQUE 
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Nose blunting is accomplished with a two point computation. The assumption 
is made that a hemispherical geometry is always maintained during ablation, and that 
the ablating material is uniform throughout. The heat entering the stagnation point is 

ae found by means of the Kemp stagnation heating Eq. (32). 
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When the temperature of the outermost segment reaches ablation temperature, 
a new radius is calculated by the following equations using an effective heat of ablation 
to model surface recession. The effective heat of ablation formulation allows for a 
reduction of that parameter as a function of pressure or shear. 


3.6.1 Laminar Flow 


For laminar flow, the heating to the tangent point is computed with the 


correlation suggested by Stetson (Ref. 12). After determining the heat rate Qr, n)* 
from a Ref, 11 table look-up, and Hat which is the minimum heat of ablation due 
min 


to either shear or pressure degradation at both the stagnation and tangency point, 
the following equations are used: 
pie Qeph At 


for recession, and for nose radius change (46) 
VR, p (Hab stag 


@ran (Hab 


!) 

2 —Sin (6) _ sta 

R, = 8-1 * PF | To-sin@) ~ T.0- smn ©) Fa) (47) 
tan 


Or, in the case where Stetson calcuiations are not prescribed: 
FDF = 0.637 DF (48) 
and Eqs. (51) and (52) are used. 


3.6.2 Turbulent Flow 


For turbulent flow, the heating at the sonic point is compuied vsing a 
correlation suggested by Hankey (Ref. 11), which is; 
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Equations (33) and (46) are again employed, with the following modifications to 
their results to reflect geometric changes. 


DFF = 1.0 - DF/R (50) 


176001941, \26000 F = 129.5 /** 


VR >; Hyp 


p. \8 
DF E 68 ( oe | (R)°* i 
FDF = 1.0 - sl (51) 


R 


where 


(50a) 


FD = (1.0 : EB), 74 (51a) 


Then the new nose radius is computed: 


n =p . FDE +2 DEF FD+DFF ea} 
n~ Bay 2.0 (DFF - FD) 
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SECTION IV 
SUBROUTINE DESCRIPTION 


4,1 INPUT 
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This is the routine for listing the data cards and then ‘reading them in according 
to NAMELIST "&BD06GS."" Input also makes minor adjustments on those data, such 
as filling in the QCQNS table with 1.0 when left blank, 


4,2 OUTINT 


This routine prints out the integer values, which were read in, before the 
computations start to allow easy inspection of that data input. 


4.3 OUTPUT 
This routine, following OUTINT, pints the rest of the data, as input. 
4.4  TRAJ 


Linear interpolation for tabular input data. The program increments time 
which is the controlling element in the tabular-interpolations. The tabular data 
evaluated within this routine include: altitude, velocity, angle of attack, internal 
convective coefficient, inside temperature and QC@NS. 


4.5 ATMOS 


The atmosphere subroutine determines the free stream pressure, density and 
temperature, The 1959 ARDC atmosphere (Ref. 13) has been assumed in this routine, 
The internal programming uses altitude in meters for the independent parameter table. 
Pressure and temperature are determined by exponential and linear curve fits valid 
over a limited altitude range. The coefficients to those curve fits are reevaluated at 
33.5, 76.2, 143.5, 161.5, 241 and 275 Kft. Free stream density is computed using 
the perfect gas relationship, Internally the units are slugs /ft®, Figures 4-1 and 
4-2 show the free stream pressure and temperature as a function of altitude, 
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100 
Figure 4-1, Pressure vs Altitude - 1959 ARDC Atmosphere 
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Figure 4-2, Temperature vs Altitude - 1959 ARDC Atmosphere 
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4,6 PRANDTL-MEYER EXPANSION (PME) 


To enter the PRANDTLL-MEYER EXPANSION routine the algebraic sum of the 
angle of attack and cone half-angle must be negative and the configuration must he either 
flat plate, cone, or wedge. Therefore, to enter this routine, angle of attack (AOFA) 
must be input as a negative term with either the flat plate, wedge or cone options, The 
PRANDTL-MEYER equations used in this subroutine assume a y= 1,4, The equation 


used to determine the initial angle associated with the given free stream Mach number is: 


UNI= y= 140,3 tan! (0,408) ~ 57.3 tan 48 (53) 
where 
2 1/2 
B= |My- (54) 
i. 
The expansion angle of interest is UNE = UNI-A+a (55) 


Three curve fits were used to determine the expansion Mach number: 
1, For UNE < 50° 


~ 26.3 


AME= GNE- 681.77 


+ 0. 7895 + 0.0279 + UNE (56) 


2. For 50° < UNE < 102, 32° 


_ __= 280,09 -3 

AME = pup cuit ay + 0.8941 + 8.28 x 10° UNE (57) 
3. For UNE > 102, 32° 

AME = ~ 284.85 __ _ 4, 9199 + 8.68 x 107° UNE (58) 


UNE ~- 130.43 


Herein remains a problem. The free-stream Mach numbers are usually high, hence 
the third equation is usually used and, in fact, its limit of applicability is exceeded if 
unbounded. A SEPARATION stop should be placed upon this equation and a ratio limit 
of expansion pressure to upstream pressure imposed, Internally this limit is 
identified as "CONPR," 
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Once the pressure ratio is determined, pressure, temperature, density and 
velocity are determined using ideal gas isentropic expansion relationships with Mach 
number as the independent parameter. Reference viscosity and referenc~ density, 
needed for the heating calculations, are computed from enthalpy and pressure dependent 
curve fits identical to those used in the shock subroutine. 


4.7 SHOCK 


The control subroutine calls the SHOCK subroutine if the configuration input is 
either flat plate, wedge, or cone, and if the algebraic sum of the cone half-angle and 
angle of attack is positive. In this routine the local (edge of boundary layer-inviscid, 
noninteraction) flow field properties, velocity, pressure and temperature, are 
computed. For a flat plate at zero angle of attack or for a wedge and cone at an 
effective angle of zero degrees to the stream flow, the local flow field properties are 
set equal to the free stream properties. If the configuration has any net positive angle 
of incidence to the free stream, oblique shock flow field data are computed for the 
effective body streamline. The basic property data utilized in this portion of the routine 
were based upon the real gas work of M. Romig. (Ref. 15) The equations are power 
series using the hypersonic similarity parameter as the independ: at parameter, The 
correlations were performed by Schadt and Livett of General Dynamics, (Ref. 13) 

The hypersonic similarity parameter for the local pressure is the product of free 
atream Mach number and sine of the cone half~angle. The parameter for velocity 

is the product of reference Mach number and sine of the cone half~angle. Figures 4-3 
and 4-4 present pressure ratio and velocity difference across a shock developed from 

a wedge configuration, Figures 4-6 and 4-7 show similar results for the cone con- 
figuration, The correlation equations actually used within the program, and their range of 
applicability, are illustrated on each figure. Note that the local property data presented 
by these figures yield values at the edge of the boundary layer and not immediately 
behind the shock, Observe also that these data are for sharp bodies only or at that 
axial station where sharp cone conditions exist. This distance can be established fror 
the entropy swallowing work done by Jacobs (Ref. 14). 
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4.8 THERMO 


The THERMO subroutine computes enthalpy for the sharp cone and sharp wedge 
configurations, The enthalpy change across the shock is determined from correlations 
of M. Romig's work (Ref. 15). Figures 4-5 and 4~8 present the enthalpy data for the 
wedge and cone respectively. For blunted cones, the THERMO routine is also used to 
compute local properties (P, V, I) based upon a pressure distribution presented by 
Griffith-Lewis and isentropic expansion from the stagnation point to that prescribed 
local pressure, The isentropic expansion assumes y= 1,40. Based upon entropy 
swallowing work by Jacobs, Ref, 13, this should be changed to y < 1.40. In any case, 
the equations were programmed to determine heat transfer and estimate flow field 
properties in the overexpansion region where sharp cone heat transfer is too high. 
Within this subroutine, recovery enthalpy and reference conditions are computed for 
use in Eckert's reference enthalpy technique for high speed heat transfer. The 
computations for recovery factor assume a Prandtl number (PR) of 0.76. Recovery 
factor is evaluated as Prandtl number to the 1/2 and 1/3 power for laminar and 
turbulent flow respectively. Reference enthalpy, density and viscosity are computed 
from Eq, (4) and the following relationships: 


for 
* < 1300 
ut = 0.426 x 107? pe 498 (59a) 
p = 0,576 x10 1% 849 p. (60a) 
for 
Pr > 1300 
ut = 0,284 x 100 04228 (59b) 
p* = 0.865 x 107 fe “+584 p. (60b) 


The superscript * indicates reference conditions and the subscript 6 is used to dencte 
edge of boundary layer pror::~ty, Viscosity, enthalpy, pressure and density have the 
usual p, I, P and p symbols. 
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4.9 DIST 


This routine historically modified the input wetted distance to allow for a 
combined run of laminar and turbulent flow. The input wetted distance was modified 
to match boundary layer velocity thickness at the Transition Reynolds number. With 
the computed laminar boundary layer thickness set as the turbulent thickness at 
transition, an imaginary turbulent wetted length was computed and used in the heat 
transfer calculations, Experience at Aerospace suggested that the matching of a 
boundary layer velocity thickness did not yield proper turbulent heat transfer rates, 
It has been suggested that a forced matching of the momentum thickness would produce 
more meaningful results. However, at this moment, no such computation is performed 
within ASAH, The input wetted length is maintained unchanged, The routine is essentially 
a dummy routine awaiting a reformulation of the mixed wetted length problem. For 
spheres and cylinders, however, this routine computes their diameter. 


4.10 KEMP 


This routine is entered whenever ICONF = 3, 4 or 5 (cone, sphere, or cylinder), 
The stagnation point calculation is performed for sphere and cylinder when ICONF = 4 
or 5 respectively, When ICONF = 3 and NOSB= 1 (blunted cone), a stagnation heat 
transfer is generated for the stagnation point transient conduction problem with nose 
blunting. For the sphere, the Kemp~Riddel equation is used to compute hot wall heat 
transfer for high Mach numbers, whereas, the Sibulkin equation is utilized at low 
supersonic speeds (Mach number less than 3.0), Equations (33) and (34) compute the 
corresponding Q for spheres, 


Note that internally the machine program computes heating on a per hour basis, 
The printout however shows the heat transfer rate in the more useful per second basis, 


Eggers, Hansen and Cunningham (Ref, 16) provide the following relationship 
for the convective heat transfer to an unswept cylinder 


Qay1 = 09319 Q (61) 


sphere 


Lees (Ref. 18) suggests (2y7l/ as the proportionality constant. The Eggers 
value is used to obtain conservative heat rates. 


4-13 . 


The following equations are employed within the ASAH program when 
considering a cylinder at a prescribed sweep angle geometry (A) and angle of 
attack (a) (useful when modeling the blunted leading edge of a supersonic aircraft 
for example) 


ANG = [ cos (a)] [ sin (A) (62) 
-1____ ANG 
8= Ta (63) 
V1 -ANG" 
Qevt = 0,981 Q, phere (©°8 s)°/2 (64) 


4,11 ECKERT 


This subroutine computes convective heat transfer for the flat plate, cone, 
wedge, and blunted cone configurations, Sphere and cylinder configuration flags shift 
the calculation to subroutine KEMP, 


For the blunted cone configuration, the program computes laminar heat transfer 
convective coefficient based upon the Griffith-Lewis corr: ‘ation (Ref. 9). This 
correlation uses the nondimensionalized wetted distance as the entry parameter 
and interpolates from a heating table carried in the program, [Eq. (18)] 


92 cos @ X 
0.802 R 


If Z > 10 the program assumes bluntness effects have been dissipated by distance, 
and the program reverts to sharp cone computations, 


Under turbulent heating conditions, the blunted cone heating is computed using 
normal shock-isentropic expansion property data generated in the THERMO subroutine, 
followed by the Schultz-Grunow turbulent skin friction equation, 


For both the sharp cone and sharp wedge configurations this ‘subroutine sequen- 
tially computes the following array of terms: skin friction, Stanton number, convective 
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heat transfer based upon the Blasius and Schultz-Grunow relationship for laminar and 
turbulent heating respectively, Implicit within this technique to determine heat transfer 
is the assumption that a modified Reynolds analogy is valid for both laminar and turbulent 
flow conditions, and that for high speed flow, Eckert's reference enthalpy parameter 

still correlates the skin friction data. Under extreme flight velocities, the reference 
enthalpy becomes a more significant measure (relative to reference temperature) of the 
energy within a flow stream as dissociation and ionization phenomenon affect the stream 
temperature, The following array of equations reflects this modeling to determine 
convective heat transfer. [Eq. (2a)] 


Re* = p*U, X'/ut 

For laminar flow conditions [Eq. (2)] 
Cf = 0.664/(Re*)!/2 

For turbulent flow conditions [Eq. (12)] 


0. 370 


Cf = [en Re*]™ 084 


Assuming modified Reynolds analogy [Eq. (11) ] 
St = = (pr) "2/3 
Computing the convective coefficient and then the convection heating: 


Hy = St Us P, (g) (3600) (65) 
Qconv = Hy (i, ~ FY) (66) 
Once again the superscript * refers to reference enthalpy conditions and Re, 
Cf, St, Hy» and Q cony 2re Reynolds number, friction coefficient, Stanton number, 
convective coefficient based on enthalpy, and convection heating respectively, Note 


that the internal units for convective heating are Btu/ft”-hr, 
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For the cone heating, both the film coefficient and the heat transfer parameter 
are multiplied by a transformation factor (from Eckert, Ref. 5) to allow for the thinning 
of the boundary layer as the cone surface area increases as a function of axial length. 
These transformation terms are V3 and (2)? 234 
respectively. 


for laminar and turbulent flow 


4,12 HANKEY-SONIC POINT SUBROUTINE 


The HANKEY subroutine requires an input of ICHK = 1 and radius (RADIUS) of 
the nose, Calculations at the sonic point are completed for either laminar or turbulent 
flow. Transition criteria within this routine are based upon either an input transition 
altitude or a Reynolds number based upon free stream properties and nose radius, For 
laminar heating the following equation is employed 


1/2, . 3,15 
: p Un IL-1, \ /,\W/2 
Qsonic = (10, 800) (22) (=: 25) (, - 129.5) \R (67a) 


For turbulent heating at the sonic point 


: b, 0.87 y \3 Iy\ /4\002 2 
Qsonic ~ 5° be 1000 vk R (7D) 


The heating rate that is provided as a boundary condition for the transient 
conditiou solution is based upon an hourly rate, 


= 3600 Q 


Qconv sonic (68) 


4.13 GETCPK 


The GETCPK subroutine provides an interpolation subroutine to determine 
specific heat and conductivity as a function of temperature, The flag needed to direct 
the calculation is KFUNT. The following code delineates the options on specific heat 
and conductivity. 


= = po, ae ae See ari oh SEEPS SE Ss aT i+ al 
FS ee aE Bi EET ED a eg Ee a SS Ee a EE Se a a aa 
ADE the MSE SEH E SS fp AS eS ee ee ee ST. PSS eS - = 3: 


perenne senreernmereres wren, een ern CORNED ONE TEE REL TESTERS, 


A ae Ob ae aie 


KFUNT = 1 specific heat and conductivity are both functions 
of temperature, 


Seoiaiatista te sali aback et 


KFUNT = 2 conductivity is a function of temperature; specific heat 
is constant. 


Sed cla hE 


KFUNT = 3 specific heat is a function of temperature; 
conductivity is constant, 


KFUNT = 4 both specific heat and conductivity are constants, 


5 <KFUNT < 10 tabular input of specific heat, conductivity and density 
versus temperature, KFUNT in this case !3 the 
nuraber of points describing the tempexaiure curve, with 
an equal number describing p, Cp, K 


i Wt ite Ral el art i ten tnt ahaa Sa 


In general (1< KFUNT < 4) the specific heat and conductivity are written as 
power series expansions of temperature. The general form for the equation is 


AES REPT Lo 


| 
B C. orK = C,+C,T+C,1T7+¢,7° (69) i. 
P 0 1 2 3 | : 
with Co» Cy Cos and Cy being input constants, The values of these constants are ee 


assumed zero if omitted in the input, The specific heat and conductivity (and density 
a if KFUNT 2 5) for a given time increment are computed for the nodal temperature 
distribution of the previous time interval, 
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4,14 ENERGY 


This subroutine performs the energy balance at the flow field/material inter- 
face, It accounts for high speed convection to the vehicle surface, corrects for current 
vehicle surface temperature, and for energy radiated away to a zero degree space 

environment [Eq. (40)]. Internal conduction is evaluated in sub~outine ONED. 


. i Ty \4 
QNET =Cq Q, -|o. 174) €| F55 


2 = @conv 
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The symbol Cq may be interpreted as an adjustment parameter that modifies 
the internally computed convective heating. For example, under turbulent flow conic 


/5 


from the Mangler Conic three-dimensional transformation. If one desires to reproduce 


heating conditions, some experts do not utilize the (ay! correction factor resulting 


4 1 
this effect within ASAH, a Cq= ae may be used to cancel the programmed 


turbulent transformation factor for the cone configuration. Another example would be 
to use Cq asa factor of safety on the predicted heat transfer rates, Factors ranging 
from 1.2 to 1.5 have been used in the conservative design of reentry heatshields. 


4.15 CONSTQ 


If NPQ(NPGAS internally) is set equal to a positive number, the CONSTQ su-~ 
routine is cailed from the Main Program. This subroutine provides a heating input 
based upon a tabular input, and the program bypasses its normal convective heating 
calculations, If NPQ is a number greater than 2, a linear interpolation of the QI table 
is completed and set into the convective heating program slot, If NPQ is less than 2, 
the convective heating term is set equal to QSOLAR, an input constant, The program 
next computes net heating by subtracting the radiation term, based upon current wall 
temperature, from the prescribed input cold-wall heating component, This option is 
useful if the pregram is to be used to solve the Fourier conduction equation only, 
bypassing the norzaai 2erodynamic heating routines, 


4.16 ONED 


This particular subroutine solves the transient one-dimensional conduction 
problem using the finite-difference form of Fourier's conduction equation, Utilizing 
subroutine GETCPK, density, conductivity and specific heat for a given point may be 
obtained as a function of current nodal temperature. The subroutine is programmed 
such that heat energy may be removed from the back face of the insulative system by 
an input convective coefficient. Within the system, conduction is the only mode of heat 
transfer; internal radiation and heat generation are not allowed. 
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For a thin-skin (one nodal point), the following equation is used to predict 
temperature change, 


TINSD - T 
‘ at__\/qnet J 
T yaa = ( Pa, Cp i (% eh (HINSD) ( 1D ) + Ty (70) 


The net energy into the system (QNET) is computed in the ENERGY subroutine, 
The heat energy removed from the system is computed from an input convective coef- 
ficient (HINSD) and an input heat sink temperature (TINSD). Density, specific heat, 
_ nodal thickness and temperature are represented by p, Cp, a; and T respectively. 
"The subscript J refers to conditions at the beginning of time interval At, whereas the 
subscript J+1 refers to the temperature at the end of the time interval, 


For conducticn problems that are divided into more than one nodal point, the 
following array of difference equations are employed, 


Aerodynamic nodal point: 


- az|QNET , Tha 7 Th 


ve (J+1) 12 A evaluated at J (40) 
0. 5| 2 + et 
n n+1 
Interior nodal point: 
T +T T -T 
is n-1 n n+1 n 
Tyger = Ay a aa a1 al (| ~ md (39) 
0.5 74 + KE 0.5 K +e 
n~-1 n nt+1 n 


Internal exposed nodal point: 


TT - 
- a-1 “n TINSD - TN 
Tys+1) = A, + HINSD |e 


a a 12 i Th J (71) 
n=-1 n ; 
0.5 kK .*r 

n~i n 

The above equations were written around the nth element and were centrally 
differenced in the interior points, In this case, the subscripts J and J+1 refer to 
the temperature at the beginning and ending of the calculation time interval (At). 
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The other symbols have the same meaning as previously, with the addition of K for 
conductivity. Al represents the heat storage term for the nth element and is defined 
helow, 


At 

Al = at | (41) 
25.0 Py an Cp,, ; ‘ 

For the situation where convection to an internal sink is specified (HINSD and TINSD), 

the internal heat flux is integrated numerically so that it can be printed out at a later time. 


4.17 ABLATE 


if an ablation temperature is provided in the input, the contro! routine calls the 
ablation subroutine. This routine uses a simple effective heat of ablation “‘hermal 
model to estimate surface recession for an ablative heatshield system, If net heat 
flux into the syste:n is positive, and if the tem»erature of the aerodynamically exposed 
nodal point is equal to the prescribed ablation temperature, this subroutine computes 
an effective heat of ablation [Eq. (46)] 


Hap Hy + 7 (ly, ~ 4y) (45) ) 


The H au and 7 used in the above equation are input parameters, usually 
determined from ground test correlations, 


For common ablative matevials, two values for » are accepted, ETAL and 
ETAT, the blowing corrections for laminar and turbulent flow respectively. Two 


Bebe a moti ede at i He hs 


equations are available to determine mass loss rate. If the input flag NQI= 0, the 
ablation rate is determined by subtracting the heat energy conducted inward from 


paving 


the net heat entering the system, and dividing the remaining energy by the materials 
effective heat of ablation, With NQI > 0, the vonvective heat transfer, rather than 
net heat transfer, is used in the mass loss computation, Results are printed in terms 


Sai Me i Wt 


of WABDOT and WAB, the mass Icss per unit area per calculation time interval and 
the integrated mass loss per unit area respectively, 


For materials that are sensitive to stagnation pressure, the user has the option 
of degrading heat of ablation at a threshold value of pressure, At that pressure (input 
as PSTMAX) the program bases heat of ablation on a second set of ablation parameters 
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(HVPSTM, ETALP, ETATP), usually reflecting a material's inferior performance i 


under a high stagnation pressure environment, 


4,18 ABLUNT 


The ABLUMT subroutine option is used to predict nose radius change for a 
sphere/cone configuration. The subroutine is entered by assigning NOSB= 1, and 
specifying a radius (RADIUS) and an ablation temperature (TWALL) in the input array. 
The blunting model assumes spherical geometry in the nose region, and uses an 


os, s 
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effective heat of ablation to predict surface recession, Laminar or turbulent heating 
is used at either the sonic point or the sphere/cone tangency point, Detailed thermo- : 
dynamic computaticns are performed at the stagnation point and engineering correlations 

are employed to estimate heating at either the sonic or tangency position, The sonic 


calculation is performed if ISTETS = 0; tangency computations are completed if 
ISTETS = 1, 
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Transition may be controlled by specification of a transition Reynolds number 
or by stipulation of a transition altitude or McCauley criteria (Ref, 2) if neither is 
specified, With respect to transition, Reynolds number is evaluated at the nose sonic 
point using sonic point flow field properties (valid only when freestream Mach number 
is hypersonic), and current nose radius, The flow field properties are evaluated as 


follows: 
Ps 2s P snie = 1160 P. (72) i 
Us = Usonio = 0-801 U,, (73) ‘ 
= es ( 7 
p Pe ite 7.24 P, (74) : 
-10 , : 
Hs = Usonic = WS1Lxl0™ U,, (75) 
i 4 = (I. -U2_. )/5.01x 104 (76) ‘ 
) sonic s ~ Usonic) /> 


The Kemp~Ridell equation, discussed earlier, is used to evaluate stagnation 
heating. Sonic point heating is evaluated using a (cos 9)3/ 2 relationship for laminar 
conditions, and Hankey's equation for turbulent conditions, These relationships are 


Rn 


“ 
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| reproduced below 
LAMINAR Q = 0.687 Q, (77) 


(assumes y= 1.2) 


p 8 U 0.3 I 
. 5.0 (7) 2) Ww 
TURBULENT Q= —> #2.) (;%e) ( = tr (78) 
(Ry”* 2 ( Ps) 1000 I, 


For laminar conditions, using the Stetson option, the heating is evaluated at the 
sphere/cone tangency point, Figure 4-9, due to Stetson, shows the tangent to stagnation 
heating ratio for laminar flow conditions, For turbulent conditions, independent of 
option, computations are always performed at the sonic point, 


It should be noted that the conduction problem is solved at only one station 


‘sat 


’ (stagnation point) when the blunting option is exercised, The program tacitly assumes 
that when the ablatic . temperature is attained at the stagnation point, it is simultaneously 
achieved at either the sonic or tangent points, 


R 
= 


Since it has become evident that material thermal performance is degraded at 
high pressure and/or high shear, this subrcutine has been extended to compute effective 
heat of ablation in a manner that admits variabie material performance as a function of 
external environment, The following equation defines effective heat of ablation in the 
nose blunting subroutine, using AI as (I, ~ I) 


Happ 2%) (FV) [HY + n, as + (%) (1 FV) [ur + ny ar| 


+ (2+ X,) (FV,) [ av, +n, a] (79) 


The subscripts 1 and 2 refer to the constituents of a two component material 
(for example the silica and phenolic of OTWR). The weight fraction of component one 
is designated Bye FV represents the fraction of material vapcrized and is input as 
a tabular function of pressure and shear. The program user may adjust material 
performance by operating on the fraction vaporized term. The terms, HV, HF, 
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Figure 4-9, STETSON - Tangency to Cone 
Laminar Heating Ratio 
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and 7 refer to heat of vaporization, heat of fusion and blowing factor. If FV and 
xX) are not input, the program internally assumes FV=1 and X, = 1, When this 


is the case, the heat of ablation equation in the nose blunting option reverts to the 
simplified form used in the ablation subroutine Eq. (45). 


When using the nose blunting option, the printout has the following meaning. 
Sonic point pressure, velocity and enthalpy are termed PSHK, VSHK and SHKI 
respectively, The stagnation heating, stagnation pressure and temperature distri- 
bution are printed in their normal array, The Reynolds number displayed is based 


upon sonic point properties and current nose radius. The stagnation point heat of 
ablation is termed HOABL. 


Additionally, the changing radius and stagnation recession are called RADIUS 
and F respectively, Further the heat rate at the sonic point (tangent point when 
"Stetson" option is used), the shear (sonic or tangent points) and the minimum heat 


of ablation (at sonic or tangent point) are also printed out and are called QSONIC, 
SHEAR and HOABLM respectively. 


With regard to degradation, the program will independently compute a J 
degraded heat of ablation for both pressure and shear conditions, and will use the 
E smaller of the calculated values to determine surface recession, 


4,19 PRINT 


q The PRINT subrou‘ine, called for all normal printing, defines those internal 
q quantities that will be printed out at those times specified by the WRITE input variable. 
4 This subroutine has no engineering interest except that se veral closed-end calcu- 


lations are performed at this juncture, Free-stream dynamic pressure, shear and 


stagnation pressure cve computed at this point. The following equations are utilized. 


; steal 
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DYNP = p, US, /2 (80) 
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QCONS is an input multiplicative factor esed to adjust shear for various 
reasons (proturberances for example). 
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Stagnation pressure is computed assuming Y= 1,4. Units for stagnation 
pressure are atmospheres, 


P 1.20 M2 |—2— | Pi, /2116 (82 
= ( bd ) ) 
STAG My (7) (2. 1) aad 


The variables printed are: TIME, ALTITUDE, U,,, angle of attack (a), Q; Hi. QNET, 
Reynolds number, time step, A weight ablated per time step, total weight ablated, 

Ts I Q asi de’ transition Reynolds number, heat of ablation, temperature of ablation, 
Ps » dynamic pressure, Us Pa | D I*, total Q, Shear, P. (atmosphexes), P. (PSF), 
emissivity, and A weight ablated per second. The calculation interval number is 
printed on the tear strip. 


In addition, if the nose blunting option is used, the radius (in feet and inches), 
estimated nose recession, Qsonic? shear at nose, and heat of ablation at surface are 
printed. 


All this is followed by a list of temperatures at each node, 
4.20 GECPK2 


This is a two-dimensional version of GETCPK, for finding Cp, K, p. Recent 
modifications now enable the engineer, as an option, to input these as tabular functions 
of temperature, as described under GETCPK, 


4,21 TWOD 


The two-dimensional conduction routine, similar to ONED, is discussed 
under 2.4.2, 


4,22 PRINTH 
PRINTH is used to print the Hollerith name of the type of configuration used. 
4,23 FINTRP 


This is an interpolation routine, used for all linear interpolation except that 
for finding ALT, VINF, AOFA, HINSD, and TINSD in subroutine TRAJ. 


4,24 AERSOP 


Prints the large Aerospace Corporation symbol at the beginning of each run, 
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SECTION V 
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INPUT DESCRIPTION 


Input to this program is accomplished by use of IBM supplied utility called 
NAMELIST (Ref. 17). The data are read from input cards, in subroutine INPUT, 
and placed on an auxiliary storage device to enuble the listing of separate data cases. 


Below (Section 5.1) is a description of the method of entry for such items 
as; Title Cards, Variables, and End of Data cards, 


Following this (Section 5.2) is an alphabetic list of variables used in one- 
dimensional problems, followed by those used for two-dimensional problems. 
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Section 5.3 then shows typical input examples so that the reader may see 
how all these tie together. 


5.1 DESCRIPTION OF METHOD OF INPUT 
i) 5.1.1 Title Card(s) 


See "TITLE" in one-dimensicnal inputs section. (5,2.1) 


<8, 


5.1.2 Variab?e Entry 
5.1.2.1 One-Dimensional 


aR at lk 


To begin data entry, "&BDO06@S" is entered in columns 2 - 8, followed by 
at least one blank, The variables and their values then follow, with "="' between the 
variable and its corresponding values, and a comma (, ) following each value, The 
values corresponding to variable names beginning with any of the letters, I, J, K, L, 
M, or N, are not entered with a decimal point although one is not necessary for other 
numbers which have no fractional part. If two or more consecutive values are the 
same, enter a fixed point number telling how many, followed by an asterisk (*) and 
the value. Following is an example of entry: 


EXAMPLE: Al = 30., 25., 25., 25., 25., 30.7 
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The data card would look like this: ‘ 
Al = 30, 4*25, 30.7, or Al = 30., Al(2) = 4*25, A1(6) = 30.7, ie 


NOTE: When entering two-dimensional arrays, such as CPTBL is 
(if the element numbers are not specified) the first element ? 
number (subscript) varies the most rapidly, i.e., for H 
CPTBL they would vary as: (1,1), (2,1), (3,1)...(1,1), : 
(1,2), (2,2)... (I, 2)......(I,d) ; 


It would be well to note, also at this point, that once a value is entered for a given 
variable, that value will remain in the computer memory, until changed. Should the : 
same value be entered more than once within one data case, however, the last value = 
entered will be used. This particular feature makes this program particularly 
applicable to running parametric studies. 


A "SEND" entered in columns 2-5 terminates the data case. 


5.1.2.2 Two-Dimensional 


For two-dimensional problems the AA, BB, CP2, DK, T2, NAR@W, and i : 
NBCGL values follow a "&TW@DIM" entered beginning in column 2, As above, equal fea 
signs and commas separate the variables and values and a "&END" in columns 2-5 
terminates each data case. For multiple cases a "&BDO6@QS" card and a "&END" 
must precede each "&TW@GDIM" card, 


5.1.2.3 End of Data 


To end each set of cases place, on the last card, "ENDDATA" beginning 
in column 2 with no data following it. Should a comment follow, this comment must 
not begin before column 9, The card images for each set of cases will be printed 
preceding the processing of each case. To terminate all processing, place two 
"ENDDATA" cards together. For further clarification, see Section 5,3 of this 
writeup and refer to the bottom of Aerospace Form 3703, "4 FIELD NAMELIST INPUT." 
For easiest data entry for the novice, Aerospace Form 2944, "80 Column FREE 
Keypunch Form" may be used. However, for orderly results put variable names in 
columns 2, 22, 42, or 62 as indicated on Form 3703, 


j i 
5-2 i) 


5.2 LIST OF INPUT VARIABLES (ALPHABETIC) 


5.2.1 One-dimensional Problems (Begin with ''& BD06QS") 
The equation for computing the conductivity [ Btu/(ft °F hr)] for node J, is 


as follows: 


K, = Aly +A2,T + A3,T" +A4_.T° where T is temperature [°R] 


or entered separately as: 


J 
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Al(J ‘ Corresponds to Al J above 
A2(J) Corresponds to A2 J above 
A3(J) Corresponds to A3 J above 
A4(J) Corresponds to A4 J above 
ALT(N)” Altitude [ft] corresponding to TIME(N) 


ALTCM(NT)®? Altitude (cm] defining a new atmosphere (see TFB) 
i ANGLE Cone or wedge half angle (degrees ] 


AGFA(N) Angle of attack (degrees] corresponding to TIME(N) 


ee ee eee ed 


The equation for computing the specific heat [Btu/lb °F] for node J; the equation ‘ 


used is: 
CP, = Bl. + B2.T + B3.T* + B41? 
J J J J J : 
where T is temperature (°R) F 


iow Palais 


or entered separately as: 


Bl(J) Corresponds to Bl J above 
B2(J) Corresponds to B2 J above 
B3(J) Corresponds to B3 J above 
BA(J) Corresponds to B4 J above 


lig ) (is less than or equal NSEG) refers to nodal distance dependent variables. 


2N) (cannot be greater than NTRAJ) refers to time dependent variables. 
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SinT) {cannot be greater than TRVRSL) refers to altitude dependent atmospheric 
=. variables. 


J and N will be used in this same context in the rest of this writeup. 
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Cu 
CONPR Limiting cone pressure ratio for Prandtl Meyer 
Expansion, The program will not lower the pressure 
ratio below this constraint. 
CPTBL(I, a4 Specific heat of segment J corresponding to TEMTBL(I, J) 
{Btu/Ib °F] 
DELTA Time step size [sec]. This is no longer a fixed increment. 
After program start-up the maximum stable increment 
will be used. A number must be entered. DELTA will 
not exceed WRITE, however. 
EKTBL(I, x4 Thermal conductivity of segment J corresponding to 
TEMTBL(I, J) [Btu/(ft °F hr)} 
EL(J) Thickness [inches] of the Jth segment 
EMIS Surface emissivity 
EMTA B(L)® Emissivity table, each value corresponding to 
TWTAB(L) ) 
ETAL Eta laminar and Eta turbulent, used in the calculation 
ETAT of the heat of ablation (see HV), a blowing correction 
ETAL2 Same as ETAL or ETAT, but for the second component 
ETAT2 within a composite (see HV2) 
ETALP Same as ETAL and ETAT, but used only when PSTAG 
ETATP is greater than PSTMAX 


4(1yMax. of KFUNT, minimum of 5, refers to individual TEMTBL vs RHOTBL, 
EKTBL, CPTBL points for material properties interpolation. 


NOTE: Only the properties of the first segment of a given material need be 
entered for the TEMTBL vs RHOTBL, EKTBL, CPTBL option, provided 
all were initialized to zero, which they are for the first data case, 


5(L) maximum of 20, refers to individual TMTAB vs EMTAB values (wall temperature 
°R vs emissivity), when emissivity is variable. (See NEMIS) i 
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FV1S(K)° Fraction vaporized component 1 (shear dependent) of 


two component composite material (see XI). 
FV2S(K) Fraction vaporized component 2 (shear dependent) 
entered as tables (max. of 10 points) vs SHRTBL(K). 


FV1P(K) { Fraction vaporized component 1 (pressure dependent). 
FV2P(K) { Fraction vaporized component 2 (pressure dependent) 
| entered as tables (max. 10 points) vs STGTBL(K). 


HEDGE(N) Input table from Entropy Swallowing Program when 
INPCW > 1, corresponds to SHKI (I; internally) 


Heat of fusion for component 1, if not all vaporized 
(see XI) [Btu/Ib] 


HINSD1(N) _ Inside heat transfer at TIME(N) [Btu/hr ft? °R] 
H@GABL Heat of ablation [Btu/Ib]} 


HREC@V(N) Corresponds to HSTAG in Entropy Swallowing Program 
(table input vs TIME1 of INPCW > 1 and corresponds 
to RECI inside this program) 


Heat of vaporization for first component of a heatshield 
{ Btu/Ib] 
Heat of ablation for second component of a composite 
heatshield 

HVPSTM HV for PSTAG greater than PSTMAX [Btu/Ib] 


ICQONF = 1 Flat plate calculation 
= 2 Wedge 
= 3 Cone or blunt cone (see NGSB) 
= 4 Sphere (see NG@SB) 
= 5 Cylinder 


5 (K) Max. of 10 refers to pressure dependent variables. 


THANK 


INPCW 


ISTETS 


NDBUG 


Hankey sonic point flag if IC@NF = 4 
Set to 1 to use, set to 0 to block usage. 
RADIUS must also be specified. 


If greater than zero, signals that cold wall inputs 
(TIME, QINT, HEDGE, UEDGE, PEDGE, and 
HRECG@QV) will be processed and used for computing 
QSUBI instead of Aero-Heating. ‘The value of 
INPCW tells how many points in the above tables 
(maximum of 50). 


For Stet on type laminar nose blunting heating 
when positive. 


Case identification number (-99 < KASE < 999) 


Determines if K and/or Cy are functions of 
temperature 


= 1 both are functions of temperature _) 
= 2 K only is a function of temperature 
=3 cy only is a function of temperature 


= 4 neither is a function of temper ‘ve, but 
both are constant. (see Ala: Bl) 


> 5 is the number of points (no more than 10) inp, 
K, C p vs temperature tables, from which 
material properties are interpolated, for each 
segment at their particular segment temperatures 
(see TEMTBL, RH@TBL, EKTBL, CPTBL) 


Debugging flag, gives a print each time through one of 
the routines listed for certain sections of each 


= 1 PRINT, ABLATE, GETCPK when 5 < KFUNT < 10 


>1ABLUNT, ABLATE, INPUT, KEMP, PRINT 
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NEMIS 


NGPSTG 


N@PSHR 


N@SB 


NPQ 


>1, #3 TRAT 

= 2 ABLUNT 

> 2 DIST, GECPK2 

= 3 PMF, ATMGS 

= 4 SHOCK 

= 5 THERMG 

= 8 HANKEY 

= 14 ENERGY/ENCO 

> 14 TW@D 

= 1 for 1-dimensional problems 
= 2 for 2-dimensional problems 


Normally zero, used as number of entries in 
variable emissivity table (EMTAB vs TWTAB), 
maximum = 20 


No. of entries in FV1P and FV2P vs STGTBL tables, 
maximum = 10 


No. of entries in FV1S and FV2S vs SHRTBL tables, 
maximum = 10 


Used with ICGNF = 3 or 4, Set = 1 for Nose Blunting. 
RADIUS and TWALL must be entered. Set = 0 if no 
blunting. 


Set to one of the positive numbers below for NON- 
AERO Heating Option 


= 1 QSUBI = QS@LAR (constant value that must be 
input. ) 


= 2< NPQ < 100, QSUBI will be interpolated from the 
QI vs TIME tables from NPQ (number of) values. 
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NPQS 


NPU 


NQI 


NSEG 


NTRAJ 


PitDGE(N) 


PFB(NT) 


Normally zero, used as number of entries in 
QSTAR vs PSTAGN tables (10 maximum) 


Punch flag - Set positive to punch cards that 


CR ie 


serve as input for Char Ablator Program. These 


punched inputs are: 


TIME = TST@R 
QSUBI = Q¢ 
PS = P@ 
RECI = HSTAB 
SHKI = HLTAB 
VSHK = UETAB 


Since a value is stored for each of the above variables 
at time of print, it is recommended that 


STOP 
WRITE > 50-NSEG to be ablated 


so that NTMP, NTME, N@P, NHS, NUE, NHL@C 
(No. table values) will be less than 50. 
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When set positive will cause the ablation routine to use 
QSUBI in place of QNET in the weight ablated (WABDG@T) 
calculation. If NQI = 0, QNET will be used iu the 


calculation, 
Number of segments in heated material, (1 < NSEG < 40) 


Number of points in trajectory table (1 < NTRAJ < 100) 
for TIME vs VINF, vs ALT, vs AGFA, QC@NS, 
TINSD, HINSD, etc. 


Table input from Entropy Swallowing Program when 
INPCW < 1 corresponding to PSHK [Atmosphere ] 


Pressure [dynes/em”} corresponding to TF B(NT) and 
ALTCM(N1) 


eee 


PSTAGN(K) “ Table of stagnation pressures corresponding to 
QSTAR(K) (see NPQS) 

PSTMAX Maximum value for PSTAG before switching from 
equation 46 to HGABL = HVPSTM x nx (I, ~ ly) 


QC@NS(N) Table of multiplying constants vs TIME used in the 
Aero Heating QSUBI equation. If QC@NS (NTRAJ) = 0, 
all QCQ@NS will = 1. 


QI(N) Table values, at Time = TIME(N), for QSUB! when 
NPQ > 1. {Btu/ft?-sec} 

QINSD(N) Back wall Q vs TIME(N) 

QINT(M) Cold wall Q from Entropy Swallowing Program 


(vs TIME) to be used when INPCW > 1. This is 
converted to hot wall Q internal to the program. 


QSGLAR A single value for heating, either above 295274. 7 ft 
(90, 000 meters) or if NPQ = 1 below: that altitude. 


QSTAR(K) Used as HV, interpolated from PSTAGN(K) 


RADIUS Of sphere, cylinder, or of blunt nose when NOSB 
= 1 [ft]. Otherwise use XDIST for sphere or 
cylinder diameter, With NOSB = 0, RADIUS will 
trigger the Griffith-Lewis Heating for cones. 


RHGQ(J) Density of Jth material segment {Ib/ft? } 


RHGFB(NT) Density of atmosphere {g"n/em*] at ALTCM(NT) 


RH@TBL(I,J) Material density of segment J corresponding to 
TEMTBL(L, J) {Ib/ft?} 


START Starting time of desired trajectory if other than zero, 
may be greater than TIME (1) [sec] 


1% has a maximum value of 10. 
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STGTBL(K) independent parameter in vaporization vs pressure 4 

table [atmospheres ] 3 

SWEEP Sweep angle [deg] for cylinder with axis not perpendicular ; 

to velocity vector : 

TABL Temperature of ablation (°R]; nose blunting requires § 

an additional input, TWALL i 

TEMP(J) _ Initial temperature of node "J" [°R] ; 

i TEMTBL(I, J) Temperature table values (when KFUNT > 5 >10) 4 
; from which material properties of segment J will be é 
. interpolated [°R] j 
TFB(NT) Temperature [°K] at ALTCM(NT)® (See TRVRSL) : 
TITLE A - 1 to 72 alphanumeric character title which is printed 

preceding each case. Must begin and end with an apostrophe é 

i 

: TIME(N) Time corresponding in trajectory to VINF(N), ALT(N) } g 
4 TINSD1(N) Temperature inside [°R] at TIME(N) 3 
TRANLT® Transition Altitude [ft] only if TRANR = 0 : 
TRANR Transition Reynolds number ‘ 

TRVRSL No. of points defining new atmosphere (see TFB) é 

4 TWALL Temperature of Ablation if NGSB = 1 [°R] must be i 


input each case or TARL will be used, i 


TWTAB(L) Temperature at wall corresponding to EMTAB(L) 


Sif any two consecutive TFB values are zero, the program will use standard $ 

atmosphere during that portion of the trajectory when two or more TFR values i 
are positive, an interpolation for TINF, RI@INF, PINF vs ALTCM is 
‘3 accomplished, but only between positive TFB values. 


a both TRANLT and TRANR = 0, McCauley Transition Criteria will be used. 


“—_ 
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UEDGE(N) Table input from Entropy Swallowing Program when 
INPCW > 1 corresponding to VSHK internally 


VINF(N) Free-stream velocity at TIME(N) ifps] 
WRITE Print interval for problem (sec] 


XDIST The distance back from the leading edge of a wedge or 
cone at which temperature profile is desired, or the 
diameter of sphere or cylinder {ft}. 


Fractional portion of compunent i of a two component 
coniposite ablative material (see ETATZ, ETAL2, 
HF, HV2, FV1, and FV2) XI taken to be 1, unless 
otherwise specified. 


5.2.2 Two-Dimensional Problems (Begin with "& TWGDIM") 


Following the normal one-dimensional inputs, the additional inputs needed 
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REE 


for two-dimensional problems (beginning with "& TW@DIM" instead of '&BD06GS") are: 


AA(J) Segment thickness (in Y direction) corresponding to 
EL(J) [inches] 


3 


BB(M) Segment width (in X direction) [inches] axial flow 
direction 


i bia - eh ats ae 


CC(J, M) > 1 for conduction into and out of grid point 
= 0 for no conduction, such as perfect insulators 


CP2(J) Corresponds to B1(J), specific heat capacitance !? 


DK(J) Corresponds to Al(J) conductivity 


T2(J, M) Initial temperature distribution in two-dimensional 
grid [°R] 


107,0-dimensional problems are assumed to be homogeneous in the X, (Parallel to 


surface) (Subscript M) direction, but not necessarily in the Y (in depth) (subscript J). 
J and M for two-dimensional problems may not exceed 10. 
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NAR@W Number of rows of AA (Max, = 10) 
NBCGL Number of columns of BB (Max. = 10) 
5.3 ZXAMPLE INPUTS 


Following are examples for inputting both one-dimensional (5.3.1) and 
two-dimensional (5,3.2) problems. These problems show certain features and uses 
of the program, The foliowing is a general description of each of the problems. 


CASE 1 is a one-dimensional, single segment problem using fixed material 
properties and deriving its heating fror. aerodynamic criteria which come from an 


input trajectory using McCauley transition criteria. 


CASE 2 is a one-dimensional, five segment problem, with C_ input as a 
function of temperature, while p and K are held constant. Heating, for this problem, 
comes from the cold wa!l input parameters associated with "INPCW = 1", This 
problem uses altitude for transition criteria. 


CASE 3 is a modification of Case 2, and adds another segment of different 
material to the face of the five segments, thus making it a six segment problem with . 
all the features described in the above paragraph, 


CASE 4 shown on the printed listing (as made by the ASAH Program) as 
Cases 4 and 5, is a two-dimensional! problem with segment 2 in slab 2 not present. 
The "START" feature was employed, allowing the program to start at some time 
other than at the first time entry. The second 'T2" card, it should be noted, supersedes 
the first "T2" card, The values on the second "'T2" card were obtained from the values 
on the first 'T2" card after running the problem for 68.7 seconds. All other parameters 
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d are as described in Case 2, above. 
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The following is a listing made by the ASAH program illustrating both the 
actual input cards used, and the standard output as given by Subroutine INPUT. 
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SECTION VI 
SAMPLE OUTPUTS 


Five cases are shown on the following pages. These correspond to the example 
input cards shown in Section 5.3 of this writeup. Descriptions of these cases are as 


follows: 


‘ere we wily Bene Hania 1 sh maaomtaiteabe tal GLUES odd Saatanaitnand 
chica aaa 


6.1 One Dimensional, One Segment, Aerodynamic Heating 


a ee 


6.2 One Dimensional, Five Segments, Non-Aeroheating 


6.3 One Dimensional, Five Segments, Non-Aeroheating, Q Correction 


a yh a Mi eels de® 


6.4 One Dimensional, Six Segments, Non-Aeroheating, No Q Correction 


ve aime tit tap Reread nee bt Sia deb aha at rats oat ee tl dt ba iesheke tetebete tet a Sn kbeatter dune 


6.5 Two Dimensional, Five Segments x Two Slabs, Non-Aeroheating 


Useful comments about the variables which are output, and what they mean 
may be found in Section 4.19. 
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SECTION VII 
RECOMMENDATIONS 


71 PROGRAM PHYSICAL MODIFICATIONS 


During the program checkout process it became evident that with a variable 
degree of effort, the physical modeling within ASAH could be expanded to enhance 
the flexibility and performance of the computer code, Toward this end, the following 
list of recommended program changes was compiled. Several of these recommen- 
dations have already been implemenfed within ASAH and are so indicated, 


7.1.1 Nose Blunting Option 

The nose blunting option has been expanded to include both laminar and turbulent 
considerations at either the sonic or tangency points. The computations are still 
constrained by assumed sphere/cone geometry, These modifications have been 
completed. However, the ultimate solution within this option would eliminate the 
spherical geometry restriction that currently plagues this subroutine, It is recom~ 
mended that the nose blunting option be further expanded to include calculations at 
many points along the nose configuration and to utilize polynomial curve fitting 
techniques to the ablated nose surface. 


7.1.2 Cylindrical Coordinates 

It is recommended that a cylindrical coordinate system option be included in 
the conduction subroutine (ONED) to complement the present rectangular coordinate 
system, The cylindrical coordinate system would be particularly useful in the nosetip 
region of a reentry vehicle, 


7.1.3 Surface Temperature 

The present ASAH computes temperatures at the midpoints of the nodal network. 
Since incipient ablation, reference flow field properties and the heat balance at the 
surface are dependent upon surface temperature, it is advisable to provide as accurate 
an estimate for surface temperature as possible. Therefore, it is recommended that 


“% 

! 
| 4 
Re 

a 
oy 
BY 
7 


= eta wet Ne ty Sait ai Sa letndhabahbatiay 
ne Sate e RRS TA cite wl ath ot Ve Le bo Ue 


DESERT AE eee rete a A Berd as Sed AES ate td 


iat Reha et Ee SSN te pes aie 


aed 


+ 


Pel is, 


Jock. 


4 


CAS eet BE AEaL Mb beta OES MARCO Dietary A ee be Bhan athe Xt 


PAM toes ay ts at 
ont 


Pcamlineue’ cake picatee anal deiteatinds nee ed aaicontiaitamatomntattaddanaateaa naka’ 


TE PREG MTT US Oe Me TE TE TES a PO FIAT SAY TEAL 


some extrapolation or analytical technique be incorporated into this program to 
allow determination of the surface tempereture from the nodal midpoint. 


7.1.4 Angle of Attack 


At present ASAH accounts for angle of attack effects by simply adding the 
input angle of attack to the body angle to determine an effective flow incidence angle. 
Although this additive. technique is adequate for very small angles of attack, it 
deteriorates rapidly as angle of ettack increases, It is recommended that a more 
sophisticated angle of attack routine be included, 


7.1.5 Ablation Temperature 


At present, ASAH performs ablation computations holding surface tempersture 
at a prescribed limiting value. It is recommended that the program be changed to 
accept a variable wall temperature as a boundary condition, 


7.1.6 Mixed Flow Boundary Layer Growth 


Heat transfer computations are performed assuming an all laminar or all 
turbulent boundary layer growth. Under actual flight conditions the boundary layer 
growth will be, at some point, determined by a partial run of laminar and turLulent 
flow conditions, The computer code should be modified to allow for this contingency 
by matching displacement or momentum thickness at the transition station, 


7.1.7 Prandtl Number 


The computer code performs calculations using a constant Prandtl number. 
Recent work by H. Jacobs of Aerospace has indicated that Prandtl number is a function 
of stream energy. Jacob's Prandtl number correlation should be added to ASAH. 


7.1.8 Plotting Option 


At present, the code does not have a plotting capability, Interpretation and 
assimilation of heat transfer results (heat transfer rate, mass loss rate, surface 
temperature, and backface temperature) could be improved if a plotting option were 
included in the computer code, 
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7.1.9 Insulation Thickness Selection 


The usefulness of ASAH would be increased if the controls of the program were 
expanded to include a recycle feature that sized a heatshield based upon perhaps a 
backface temperature design criterion, With this option a heatshield thickness could 
be determined with a single computer simulation, 


7.1.10 Wetted Distance Adjustment 

For computations near the nosetip, wetted distance may significantly change 
as nosetip recession proceeds to grind away at the RV body, It is recommended that 
a change be made to ASAH to provide for variable wetted distance as a function of 
nosetip recession, This particular option would be most useful in the strong entropy 
swallowing regime near the RV nose. 
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